x_length = length(L0s);
X=L0s;
Y=zeros(x_length,1);
for i=1:length(L0s)
    Y(i)=Ks(1,3,i);
end
plot(X,Y);
hold on;

for i=1:length(L0s)
    L_i = L0s(i);
    K=L2K(L_i);
    Y(i)=K(1,3);
end

scatter(X,Y);



% xy_length = length(L0s);
% [X,Y] = meshgrid(L0s,L0s);
% Z=zeros(xy_length,xy_length);
% for i=1:length(L0s)
%     for j=1:length(L0s)
%         Z(i,j)=Ks(1,1,i,j);
%     end
% end
% surf(X,Y,Z);
% hold on;
% 
% for i=1:length(L0s)
%     for j=1:length(L0s)
%         L_L = L0s(i);
%         L_R = L0s(j);
%         K=L2K(L_L,L_R);
%         Z(i,j)=K(1,1);
%     end
% end
% 
% scatter3(X, Y, Z);
% hold off;